Magnetic-Field-Induced Antiferromagnetism 
in Two-Dimensional Hubbard Model: Analysis of CeRhlns 
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We propose the mechanism for the magnetic-field-induced antiferromagnetic (AFM) state in a 
two-dimensional Hubbard model in the vicinity of the AFM quantum critical point (QCP), using 
the fluctuation-exchange (FLEX) approximation by taking the Zeeman energy due to the magnetic 
field B into account. In the vicinity of the QCP, we find that the AFM correlation perpendicular to 
B is enhanced, whereas that parallel to B is reduced. This fact means that the finite magnetic field 
increases Tn, with the AFM order perpendicular to B. The increment in Tn can be understood in 
terms of the reduction of both quantum and thermal fiuctuations due to the magnetic field, which 
is caused by the self-energy effect within the FLEX approximation. The present study naturally 
explains the increment in Tn in CeRhlns under the magnetic field found recently. 
KEY WORDS: field-induced magnetism, FLEX, antiferromagnetic fluctuations, CeRhlns 



Recently, critical phenomena in the vicinity of the mag- 
netic quantum critical point (QCP) have attracted much 
interest in strongly correlated metals. Experimentally, 
the outer magnetic field is frequently used to change 
the distance from the QCP. As for the antiferromagnetic 
(AFM) QCP, the magnetic field is believed to increase 
the distance to the QCP in general. Spin fluctuation 
theories such as the SCR theory [1] and the fluctuation- 
exchange (FLEX) approximation [2], have succeeded in 
describing various critical phenomena in metals close to 
the AFM-QCP, such as the non-Fermi liquid-like behav- 
iors of various transport coefficients [3,4]. However, pre- 
vious studies on the effect of the magnetic field based on 
the spin fluctuation theory were not comprehensive [5,6]. 

CeMIus (M=Rh, Co, or Ir) is a well-known quasi-two- 
dimensional heavy fermion compound, where single con- 
ductive Cein layers stack perpendicular to the c-axis. 
CeCoIns is a superconductor with Tc = 2.3 K at am- 
bient pressure [7]. In CeRhlns, the AFM order emerges 
at Tn = 3.8 K at ambient pressure, and the supercon- 
ductivity emerges at Tc « 2K below P = 1.6 GPa [8,9]. 
Recent experiments reveals that the Tn increases under 
the magnetic fleld along the a(6)-axis. When B = 9 T, 
the increment in Tn is approximately 0.15 K. A small 
increment in Tn is also observed in Ce2RhIn8 which is 
composed of double CeIn layers. However, there has been 
no theoretical explanation for this phenomenon. 

In the present study, we investigate the two- 
dimensional Hubbard model under the uniform magnetic 
fleld B along the x-axis, based on the FLEX approxima- 
tion. In the vicinity of the AFM-QCP, we flnd that the 
AFM spin correlation of the y(z)-component is enhanced 
by the applied magnetic fleld. In the obtained phase di- 
agram, the magnetic transition temperature Tn, below 
which the staggered magnetism emerges on the j/z-plane, 
increases with magnetic field. The mechanism of the 
field-induced antiferromagnetism (FI-AFM) proposed in 
the present study will be universal in low-dimensional 
metals close to the AFM-QCP, contrary to the fact that 



the magnetic field suppresses Tn in usual models by the 
mean-field approximation. The present study naturally 
explains the enhancement in Tn under the magnetic field 
in CeRhlns. 

We analyze the following two-dimensional Hubbard 
model: 
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where a = 1(— 1) corresponds to the t- (i-) spin state and 
Eko- = Ek + crB, where the factor aB represents the Zee- 
man energy. The spin quantization axis is the a;-axis. 
We study the square lattice tight-binding model with 
nearest neighbor hopping (t) and next-nearest neighbor 
hopping {t'). The dispersion of the electron is given by 
Ek — — 2t(cosfc^ + cosfcy) — At' cos kx cos ky . We study 
the case of {t,t') — (1,-0.25) with the electron density 
n = 0.90 {n — 1.20) per site, which corresponds to a hole- 
doped (electron-doped) high-Tc cuprates. In the case of 
n = 0.90, the Fermi surface(FS) is very close to the van- 
Hove singular point (at (tt, 0) in this case; see Fig. 3), 
and it is similar to the largest (main) cylindrical FS in 
CeMIns (M=Co,Ir,Rh) [10]. Assuming a similar single 
cylindrical FS, many aspects of CeMIns, particularly the 
(ix2-y2-wave superconductivity, can be reproduced by the 
perturbation study [11,12]. 

In the presence of the magnetic field along the x-axis, 
the dynamical spin susceptibilities within the FLEX ap- 
proximation (or random-phase approximation (RPA)), 
xUq) and xl(z)iq), are given by 

Xyiq) = Xz{q) = iXTdiq) + Xuiq)) /4 (2) 

xAq) = [xt,t(9) + Xi,i('?)]/4 + Uxtdiq)xliiq)/2, (3) 
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Xl,, {q) = -TY,Ga{k + q)G^' (fc) . (6) 

fe 

Note that Xtdil) — {Xi,T('i')}*- Here and hereafter, 
we promise that q = (q, iw„) = (q, 27r«nT) and k = 
(k, ze„) — (k, 7rz(2n + 1)T). Apparently, both Xx{<l) and 
Xy(z){Q) are even functions oiB, reflecting the reflectional 
symmetry in spin space. Apparently, Xx (?) = Xy (<?) when 
B = 0. 

The self-energy in the FLEX approximation is given 

by 

E^ik) = U^TY)Ga{k - q){X-a,-a{Q) - X-a,-aiQ)) 
Q 

+G-„{k - q)Xa,-a{q)] + Un-a, (7) 

where = TX;fcImG<,(fc)e-^'=" ° / TT is the density of 
electrons with a-spin. Here, we solve the Eqs. (2)-(7) 
together with the Dyson equation G~^{k) = ie„ + /x — 
Ck — (tB — Sct(A) numerically, by adjusting the chemical 
potential so that n = n^. 

Here, we discuss the numerical results obtained by the 
FLEX approximation. We use 64 x 64 k-meshes and 1028 
Matsubara frequencies in the present numerical study by 
FLEX approximation. Figure 1 shows the obtained static 
staggered spin susceptibilities: xS"*^ = maxq (Q) 0) , 
where a — x,y,z. as = maxq C/x''(q, 0) is the Stoncr 
factor without B. In the FLEX approximation, as < 1 
is always satisfied at finite T in two-dimensional sys- 
tems, so the Marmin-Wagner-Hohcnbcrg theorem is sat- 
isfied [13,14]. The momcntimi dependence of Xa{(l,0) 
{a = X, z) and the splitting of the FS under the mag- 
netic field are given in figs. 2 and 3, respectively, in the 
case of n = 0.90. 

In Fig. 1, x™^^ decreases whereas x™^'^ increases with 
B II X in both cases of n = 0.90 and n = 1.20 by FLEX 
approximation. Their field dependence becomes more 
prominent as U increases, that is, as ag approaches unity. 
These results indicate that the distance to the AFM-QCP 
decreases owing to the uniform magnetic field. In the 
FLEX approximation, the field dependence of the suscep- 
tibility is caused by (i) the change in the nesting condi- 
tions owing to the Zeeman splitting of the FS, and (ii) the 
self-energy effect (or mode-mode coupling effect) which 
represents the reduction in y^^^ and its Curie- Weiss-like 
temperature dependence owing to the spin-fluctuations. 
In the FLEX approximation, a large ImS(k, —i5) caused 
by spin fluctuations reduces the density of states (DOS) 
at /i, which makes ;^i^lex x^^"^- Below, wc will dis- 
cuss that the effect (ii), which is absent in the RPA is 
important to explain why xZt^ is enhanced under the 
magnetic field parallel to the a;- axis. 

We discuss the physical reason for the field enhance- 
ment of the AFM correlation: First, the uniform mag- 
netization induced by B || x will reduce the AFM corre- 
lation along the x-direction. This leads to the enhance- 
ment of x™*" by contraries, as a result of solving the 



conflict between spin-fluctuations with different compo- 
nents. The increase in x™^'' will be more prominent in 
lower dimensional systems because the reduction of Tn 
due to fluctuations is large in general. Note that the re- 
duction of the staggered moment at T = owing to the 
quantum fluctuations is approximately 40%(15%) in two 
(three) dimensional S = 1/2 Heisenberg model without 
a magnetic field. 

Consistently with the above discussion, X™(2^ increases 
whereas x^'^^ decreases under B || x in the present model 
by the FLEX approximation. We have checked that this 
is a universal behavior in two-dimensional systems close 
to the AFM-QCP, by studying various types of Hubbard 
model. Here, we briefiy discuss the self-energy effect for 
susceptibilities: When B = 0, XuC^'O) by the FLEX 
approximation is reduced from the RPA's value because 
of the reduction of the DOS, which is caused by the 
large ImS under strong spin-fluctuations. Considering 
that S(fc) « Xa{q)G\k + q), the change 

in Xt,j.(Q)0) within the lowest order with respect to the 
self-energy is given by 

S'xlli^, 0) « -T^ G^kfC^ik + q)GO(fc + q') 

k,q' 

x2U^Xx{q') + 2xy{q')). (8) 

Because it is negative, x^^^ in tbe FLEX approximation 
becomes smaller than that in RPA. Once B || x is ap- 
plied, the reduction in Xx{q') owing to the field- induced 
uniform magnetization will make |(5'x| |(Q, 0)| smaller. 
Consequently, x™'^^ increases in proportion to B^ as long 
as only the self-energy effect is taken into account. 

On the other hand, unphysical results are obtained by 
RPA, where all G's in Eqs. (2)-(7) are replaced with 
Go's. In the case n = 0.90, both x™^'' and x^"""" by RPA 
increase with B as shown in Fig. 1, possibly reflecting 
the fact that the FS is close to the van-Hove singularity. 
On the contrary, both x™^'^ and x'^^^ decreases with B 
when n = 1.20. Thus, results given by the RPA are not 
universal, depending sensitively on the shape of the FS. 
As a result, the self-energy effect included in the FLEX 
approximation is indispensable in reproducing the phys- 
ically reasonable behavior of the two-dimensional nearly 
AFM metals (i.e., as ^ 0.98) under a magnetic field. 

In the next stage, we study the magnetic-field de- 
pendence of the Neel temperature Tn by assuming a 
weak three-dimensional coupling [13,14]. To simplify the 
analysis, we define Tn in the presence of the magnetic 
field under the condition that maxq [/x'J |(q, 0) = «§, 
where a§ is a constant which is slightly smaller than 
unity. By putting ag = 1 — J_\_/U ~ 0.99 (Jx denotes 
the interlayer magnetic coupling strength), we obtained 
reasonable Neel temperatures of k;-(BEDT-TTF)2X and 
TMTSF based on the dimer model [13,14]. Figure 4 
shows the field dependence of Tn given by the FLEX 
approximation, for several ag's. We find that the field 
enhancement of the Neel temperature occurs in nearly 
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AFM metals in two dimensions, which has been pointed 
out in the present work for the first time. In Fig. 4, Tn 
starts to increase in proportion to B^, and it ahnost sat- 
urates at approximately B* ~ 0.3 when n = 0.90. This 
result also means that the system approaches the AFM- 
QCP by applying a magnetic field. The increment in Tn 
is larger when n = 0.90, reflecting the closeness to the 
van- Hove singularity. 

Here, it is notable that in the antiferromagnetic 
isotropic Heisenberg chain under the magnetic field along 
thea;-axis, {SfSp-M"^ oc cos(27rM(i- 
j)) and {SfSj) oc {—iy~^\i — where rj decreases 

from unity with the magnetic field [15]. Their field de- 
pendencies are consistent with the present study of a two- 
dimensional Hubbard model. In the XXZ-Heisenberg 
chain, an infinitely small magnetic field along the a;-axis 
induces the staggered magnetization of the y-component 
in the case of < Jx [16,17]. In the opposite case, 
Jz > Jxi the staggered magnetization along the z-axis, 
which exists without the field, is enhanced by B \\ x 
[18]. We also point out that Fukusima and Kuramoto 
studied a localized electron model with interactions be- 
tween quadrupolc moments by a local approximation, 
and found the field enhancement in Tq due to the sup- 
pression of fluctuations [19]. 

Note that the field-induced SDW is realized in the 
quasi-one dimensional metal, TMTSF, owing to the or- 
bital motion of electrons, free from the Zeeman effect 
[20]. However, various characteristics of the FI-AFM in 
CeRhIn5 do not coincide with that observed in TMTSF. 
In fact, CeRhlns possesses both cylindrical and spheri- 
cal FS's. They arc naturally explained in terms of the 
Zeeman effect as discussed in the present study. 

We now discuss the experimental results of CeMIn5 in 
the present study. The bandwidth of the present model 
is ~ 10. If we estimate the renormalized quasiparticle 
bandwidth of CeMIns to be ~1000 K [21], the temper- 
ature T = 0.02 corresponds to ^2 K, which is close to 
Tc in CeCoIng [22]. The magnetic field 5 = 0.1 in the 
present work corresponds to ^ 5 T for the M = ±5/2 
Kramers doublet (KD), because the Zeeman energy for 
Ce-''+ is {6/7)fiBMH (6/7 is the 5-value of Ce^+). Note 
that the renormalization factor averaged over the FS is 
0.217 in the present FLEX approximation for !7 = 5 at 
T = 0.02. Tn in CeRhlns continues to increase with the 
magnetic field parallel to the afe-planc, at least below 9 T; 
Tn = 3.8K at OT, and Tn(9T) - Tn(OT) « 0.15 K [8,9]. 
Whereas Tn decreases monotonically when B || c, as is 
observed in usual 3D heavy Fcrmion systems. This is 
naturally understood because the orbital motion of elec- 
trons, which is absent in the present study where B is 
parallel to the 2D system, will destroy the AFM state to 
obtain the energy due to the Landau diamagnetism. 

Furthermore, we discuss the anisotropy of x{q) in 
CeRhlns: The lowest KD of Ce-'^+-ion in CeRhlng 
is rf^ \z;±) = f3\M, = ±5/2) - a\M, = t3/2) 
[23,24], which is approximately 70 K lower than the 



second lowest KD. If we put {a,f3) « (0.44,0.9) [24], 
(^;±|J^|.2;±) = ±(2.5/32 - 1.5a2) « ±1.74. On the 
other hand, {x; ±\Jx\x; ±) = ±\/5a/3 w ±0.885, where 
|.x;±) = (|z;+) ^ ))/\/2. Then, the anisotropy 
of the susceptibility of a single Ce'^+-ion is Xa/Xc ~ 
1.74/0.885 = 1.97, which is similar to the experimen- 
tal ratio. On the other hand, several neutron experi- 
ments on CeRhlns revealed that the magnetic moments 
on Ce sites lie on the afe-plane below Tn, whose effective 
moment is /ieff = 0.264/iB [23,24]. This suggests that 
the antiferromagnetic RKKY interaction between near- 
est neighbor Ce sites is XY-like; Ja.t > Jc [25]. Then, 
the magnetic field along the a-axis will enhance the AFM 
correlation along the 6-axis as a result of the reduction 
of fluctuations, which is similar to the behavior of the 
XXZ-Heisenberg chain under B [16,17]. In fact, /Ueff 
is much smaller than {Q/7)iJbB{x;+\Jx\x;+) w 0.76/iB, 
which suggests that the quantiim fluctuations arc strong 
in CcRhlns, reflecting its two-dimensionality. As a result, 
the field-enhancement in Tn observed in CeRhlns is well 
understood in terms of the reduction of spin-fluctuations 
by a magnetic field. It is a future research problem to 
take the Kondo effect into account beyond the FLEX 
approximation. 

In summary, on the basis of the FLEX approxima- 
tion, we found the field-induced antiferromagnctism in a 
two-dimensional Hubbard model, as a result of solving 
the conflict between fluctuation with different directions 
by a magnetic field. This phenomenon is expected to 
be prominent and universal in the vicinity of the AFM- 
QCP in lower dimensional systems, irrespective of the 
fact that the field-induced uniform magnetization tends 
to decrease the AFM moments. The induced AFM mo- 
ments are almost on the plane perpendicular to the ap- 
plied magnetic field, to earn the Zeeman energy by cant- 
ing. Experimentally, the field-induced increment in Tn 
will be more prominent when B is parallel to the 2D sys- 
tem, because the reduction of Tn caused by the orbital 
motion effect (Landau quantization) is absent. 

As for two-dimensional organic metals, the field- 
induced transition from the paramagnetic metal to the 
AF insulator is found in k;-(ET)2Cu[N(CN)2]C1 under a 
pressure, below Tco = 13 K and above Hc2 [26]. Also, 
field-induced SDW due to the Zeeman effect is expected 
to be realized in r-phase organic metals [27]. These phe- 
nomena will be explained by the present mechanism [28] . 

Finally, we note that the present results by the 
FLEX approximation seems reasonable in terms 
of the fluctuation-dissipation theorem; {Sa^) = 
EqC^ImXa(<z)coth(^) (a = x,y,z), and 
^x,j/,2^^^2^ ~ ^ in strongly correlated systems. When 
Xa{Q) grows, {Sa^) will increase (especially in 2D sys- 
tems). Thus, Xy(z){Q) will increase when XxiQ) de- 
creases by B II X [28]. 

We are grateful to K. Yamada, D.S. Hirashima, K. 
Miyake, Y. Kuramoto, M. Tsuchiizu and Y. Matsuda for 
useful discussions. 
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FIG. 1. Field dependences of x™''" and Xy'"'' by the FLEX 
approximation (RPA), under conditions (a) n = 0.90, U = 5 
{U = 2.10) and (b) n = 1.20, U = 8 {U = 2.96), at T = 0.02. 
In the numerical study by RPA, 256 x 256 k-meshes and 256 
Matsubara frequencies are used. 
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FIG. 3. FS's for |- and J.-electrons under various magnetic 
fields for n = 0.90 and 1.20 at T = 0.02. FIG. 4. Obtained phase diagram for Tn vs B for various 

ag's. 
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